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1.0  INTRODUCTION  j 

i 

During  the  third  quarter  of  the  Nonlinear  Dome  jj 

Effects  Study  some  numerical  results  have  been  obtained  which  j 

give  insight  into  the  influence  of  a thin  dome  upon  parametric 
generation  efficiency. 

2.0  THIN  DOME  RESULTS  j 

The  infinite  planar  dome- transducer  model  is 
shown  in  Fig.  1.  The  transducer  is  simultaneously  excited 
at  two  frequencies. 

Our  simplified  model  led  to  the  following  integral 
for  the  normalized  difference  frequency  pressure  distribution^" 


- sin  (iH-0) 

• e r#dr,  d*  . 


Fig.  2 shows  the  relationship  between  R0 , R,  r,  ro  and  the 
angles  $ ana  In  Eq.  (1)  A is  given  by 


1- 


XKL 


FIG.  2 GEOMETRY  PERTINENT  TO  Tz  AND 
THE  CALCULATION  OF  P2  («_). 
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e1  ( T + dk-  C0S  *) 


(2) 


v.v* 


i 


i 

! 


■ i 
.1 


and  z represents  the  attenuation  path  length  of  the  primary 
beams  to  the  source  point  as  shown  in  Fig.  2. 


r sin'/' 
o 

cos  e 


(3) 


In  order  to  evaluate  the  integral  in  Eq.  (1)  it 
was  first  reexpressed  in  rectangular  coordinates.  Integration 
with  respect  to  x was  done  analytically  and  integration  with 
respect  to  y was  done  numerically  using  Gaussian  quadrature. 

A representative  field  point  was  taken  at  2500  yards  along 
the  y axis  from  the  dome.  The  virtual  array  was  arbitrarily 
truncated  at  the  range 


* 

r 


1 

al+  a2 


(4) 


where  and  are  the  attenuation  coefficients  at  the  primary 
frequencies.  This  is  justified  since  we  are  only  interested 
in  relative  amplitudes  of  pressure  as  the  steering  angle  and 
dome  parameters  are  varied. 
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of  the  steering  angle  0 for  a V steel  plate  placed  10"  in 
front  of  the  transducer.  The  no  dome  case  is  obtained  by 


letting  the  dome  thickness  h be  zero  which  simplifies  and 

* 

D£  . Originally,  it  had  been  thought  that  a comparison  could 

be  made  with  the  infinite  plane  wave  solution  to  Burgers' 

equation  which  was  quoted  in  the  Semi  Annual  Report^.  This 

is  not  feasible  because  we  would  have  to  evaluate  Eq.  (1)  in 
* v 

the  limit  as  r ^<x>  . 

*2 

Table  I gives  - — computed  by  Eq . (1)  at  a point 
2500  yards  from  the  dome  along  the  y axis  as  a function  of 
steering  angle  and  primary  frequency.  Even  though  we  have 
assumed  that  the  dome  is  acoustically  thin, it  is  evident 
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Steel  Plate  Dome  -7—  Thick 

4 

11 

10  From  Transducer 

Primary 

Frequencies 

kHz 

Steering  Angle  6 

0°  20° 

45° 

50,60 

1.9  10"2 

3.5  10"5 

8.5  10"5 

40,60 

2.75  10"3 

6.5  10"6 

5.9  10‘6 

40,70 

1.21.  10“2 

3.6  10"6 

1.32  10“5 

No  Dome 

0 

it 

A 

0° 

20° 

45° 

50,60 

3.7  10_1 

1.5  10"4 

6.7  10'5 

40,60 

2.5  10-1 

2.1  10-4 

7.7  10"5 

40,70 

4.1  10"1 

2.0  10"4 

7.2  10"5 

i: 


TABLE  I.  ~ on  the  y axis  (see  Fig.  1)  at  r - 2500 


yards.  SPL  at  transducer  is  111.6  dB  ref 


1m  bar 
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A 

P 

will  vary  with  these  two  parameters. 


2. 


In  conclusion,  we  can  say  that  in  most  cases  the 

A 

p2 

amplitude  of  — = is  quite  sensitive  to  variations  in  w., 

“o 

uy , and  0 . 
x 2 


3.0 


EXTENSION  OF  THE  MODEL  TO  THE  THICK  DOME  CASE 


For  conventional  sonars  the  assumption  that  the 
dome  is  acoustically  thin  has  been  a good  one.  Parametric 
sonars  require  two  high  frequency  primaries.  Thus  the  thin 
dome  assumption  may  not  always  be  justified.  In  the  past, 
a dome  has  been  assumed  to  be  acoustically  thin  if  its  thick- 
ness is  th  of  the  radiating  wave  length.  It  would  be 
worthwhile  to  check  the  validity  of  such  a criterion. 
Consequently,  the  model  has  been  reformulated  for  a thick 
dome.  The  integral  in  Eq.  (1)  is  unchanged.  Only  the  co- 
efficient A is  modified  and  this  is  included  in  Appendix  A. 
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APPENDIX  A 

To  extend  our  plane  model  to  the  thick  dome  case  we  draw 
1 2 

upon  Brekhovskikh  and  Van  Buren  . We  retain  the  same  nomen- 

2 

clature  as  Van  Buren  . Subscripts  L and  S refer  to  longitudinal 
and  shear  waves  respectfully.  Region  1 extends  from  the  dome 
to  the  transducer .Region  2 lies  inside  the  dome  and  region  3 extends 
beyond  the  dome.  Brekovskikh^  gives  the  expressions  for  the 
reflection  and  transmission  coefficients.  We  need  to  determine 
the  amplitude  of  the  transmitted  wave  in  region  3 in  order  to 
compute  the  virtual  array.  When  media  1 and  3 are  identical 
liquids,  as  in  our  case,  then  the  reflection  coefficient  is 
given  by^ 


_R.  . 1(m2'm2-1) 

1 m + (u2  - M2  + 1 J i 


(1) 


and  the  transmission  coefficient  by 

2N 


JL 

i 


2M  + 


^N2  - M2  + 1 ^ i 


(2) 


Here  I,  R,  and  T refer  to  the  amplitudes  of  the  incident, 
reflected, and  transmitted  waves.  M and  N will  be  defined 
later.  The  boundary  condition  at  the  transducer  face  is: 
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1 

i 


0 


-ikx  8 in  9. 


Vo* 


LI 


(3) 


Here  P is  the  steady  state  pressure,  k the  wave  number, 

^ the  density,  u the  angular  frequency,  and  9^  the  angle 
of  incidence  in  region  1 upon  the  dome.  In  region  1,  we 
assume  that  P has  the  following  form: 

^ / ivy  -ii7y\  -iKx 

P - ^ Ie  + Re  J e . (4) 


Where  v « k cos  *L1  and  K - k sin  «L1  Eqs.  (3)  and  (4)  can  be 
used  to  relate  R and  I yielding 


By  Eqs.  (5)  and  (1),  I is  found  to  be 


The  quantities  N and  M are  defined  by 
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N - 

ZL2  co®2  2*S2 

. ZS2,ln22<S2 

(7) 

ZL1  ,ln  QL 

zLl«lnQs 

M - 

ft)' 

f 2 ' 

^co«  2«s2  i 

•COt  \ +(l7i)(Sln2  2#S2  , 

) cot  % 

(8) 

and  Ql  - 

d kL2 

cos  tu 

(9) 

QS  “ 

d kS2 

coa  "S2 

(10) 

where  d is  the  dome  thickness.  In  the  above  formulae  Z is  the 
acoustic  impedance  and  is  defined  by 


Z - 


oV 
cos  9 


The  velocities  for  longitudinal  and  transverse  waves  in  the 
dome  are  given  by 


and 


(Ha) 

(11b) 
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X and  m are  the  Lame  constants  given  by 


(l+')(1’2') 


(12) 


2 f l + i 


where  » is  Poisson's  ratio  and  E is  Young's  Modulus, 


By  Snell's  law,  we  have 


sin  8r 


sin  0, 


•fe) 

■(•*) 


sin  8 , 


sin  8 


LI  * 


(14a) 


(14b) 


The  wave  numbers  for  the  transverse  and  longitudinal  waves  in 


the  dome  are  given  by 


*L2  " 


(15a) 


(15b) 


The  first  digit  of  the  subscripts  denotes  the  region  (3  is 
| beyond  the  dome)  and  the  second  digit  indicates  the  primary 

frequency.  The  difference  frequency  wave  number  k.  is 
[j  defined  by 

k.  - k9  - k.  . (21) 

r 

It  is  assumed  that  w2  ^ W1  * N and  M given  by  Eqs.  (7)  and  (8) 
are  functions  of  the  frequncy  via  QL  and  Qg  given  by  Eqs.  (9) 


•ft 

and  (10).  In  order  to  find  U we  must  compute  T31  T32  from 

Eq.  (18).  From  here  on,  we  shall  omit  the  first  subscript  3 
Indicating  the  region  beyond  the  dome  and  simply  retain  the 
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second  digit  to  refer  to  the  appropriate  primary  frequency.  In 
view  of  Eqs.  (17a, b),  we  can  write 


4 A ' • " 


T1T2 


1 2 


where 


“ ‘ (ala2  + blb2)  (el*2  + flf2  ) 

- (bla2  ' Vl)  (f2el  - fle2) 


The  magnitude  of  TjTg  is 


T1T2 


2 Nx  jyj2 V”2  + ji 


9 *1 
1 2 


(22) 


(23a) 


J - ( Vl  - Vl  ) (*1*2  + blb2  ) 

+ (eie2  + fjf2  ) (v2  - b2aj  (23b) 


(24) 
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Now  we  can  define  U in  Eq.  (20)  by 


U - NjN2  j H + Ji 


(25) 


Since  we  only  want  Che  magnitude  of  A,  then  we  have  by  Eq. 


0 P„  k-  2 N2  ^H2  + J2 


'l  V1  V2 


z 


50  EVERGREEN  PLACE  / EAST  ORANGE.  NEW  JERSEY  07018 


APPENDIX  A 
REFERENCES 


1.  L.  Brekhovskikh,  Waves  in  Layered  Media 
New  York,  N.Y.,  1960),  pp.  68-69. 

2.  A.  Van  Buren  and  M.  Breazeale,  J.  Acoust 
1014-1020  (1968). 

i 


i 


(Academic  Press, 


. Soc.  Amer.  44, 


I 

* . 


17 


